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ABSTRACT 

I suggest a theoretical quantitative definition for the termination of the 
asymptotic giant branch (AGB) phase and the beginning of the post-AGB phase. 
I suggest that the transition will be taken to occur when the ratio of the dynam- 
ical time scale to the the envelope thermal time scale, Q = Tdyn/rxH-env, reaches 
its maximum value. Time average values are used for the different quantities, as 
the criterion does not refer to the short time-scale variations occurring on the 
AGB and post-AGB, e.g., thermal pulses (helium shell fiashes) and magnetic ac- 
tivity. Along the entire AGB the value of Q increases, even when the star starts 
to contract. Only when a rapid contraction starts does the value of Q start to 
decrease. This criterion captures the essence of the transition from the AGB to 
the post AGB phase, because Q is connected to the stellar effective temperature, 
reaching its maximum value at T ~ 4000 — 6000 K, it is related to the mass loss 
properties, and it reaches its maximum value when rapid contraction starts and 
envelope mass is very low. 

Subject headings: stars: AGB and post-AGB 



1. INTRODUCTION 

The transition of asymptotic giant branch (AGB) stars to the early planetary nebula 
(PN) stage is a poorly understood phase of stellar evolution. It is agreed upon that a rapid 
drop in mass loss rate accompanies the transition from the AGB to the post-AGB phase. 
However, no well defined objective criterion exists for the transition from the AGB to the 
post-AGB phase (the termination of the AGB). Blocker (1995), for example, used both mass 
loss rate and the stellar pulsation period, i.e., a dynamical property of the star. He took 
the mass loss rate to decrease with decreasing pulsating period, and defined the zero-age 
post-AGB phase when the pulsation period is 50 days. Below I will suggest a new property 
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to theoretically define the transition from the AGB to the post-AGB phase. It is based on 
both the thermal and dynamical properties of the star, and is connected to the mass loss 
rate and mass loss inhomogeneities. 

Observationally, AGB stars are well defined, and visible well developed post-AGB stars 
can be clearly defined in principle, but practically there are many difficulties (e.g., Hrivnak 

et al. 1989; Szczerba et al. 2007; Suarcz ct al. 2006). Still, there is no theoretical definition 
of the transition. Recently, Szczerba et al. (2007) classified many post-AGB stars. However, 
they classified mainly unobscured stars in the IR and/or the visible band, most likely well 
after they have left the AGB. Observationally, most stars are likely to be obscured during the 
transition because of a high mass loss rate. Several parameters were proposed to define the 
transition from the AGB to the post-AGB phase, but each one of them has some problems. 
A drop in the mass loss rate. This criterion (e.g., Suarez et al. 2006) captures the essence of 
the transition. However, what mass loss rate should be used to mark the transition? What if 
the star rotates and/or interacts with a binary companion such that mass loss rate depends 
on the companion? 

Optical depth. Suggestions have been made that the transition occurs when the totally ob- 
scured AGB star becomes visible again. But at what wavelength? Another problem with this 
definition is that it depends on the viewing angle if the mass loss geometry is not spherical. 
In addition, low mass AGB stars with low metallicity might never become totally obscured. 
Pulsation period. Although mass loss rate is tightly coupled to pulsation, it is not clear what 
pulsation period should be used. Also, AGB stars with different total masses and luminosity 
can have the same pulsation period at very different effective temperatures. The usage of 
pulsation (e.g.. Blocker 1995) includes in it the dynamical time of the star, but no 'natural' 
transition value exists. In the presently proposed criterion the dynamical time is included 
with a quantitative measure. 

A change in the dominating mass loss m,echanism. On the AGB the mass loss process 
is that of pulsations coupled with radiation pressure on dust, while for the central stars of 
PNs it is mainly radiation pressure on ions. The idea is that the transition is defined when 
the dominate mass loss process switches from pulsation and radiation pressure on dust to 
radiation pressure on ions. There are two main problems with this. Firstly, the physics 
is not well understood to connect this transition to stellar evolutionary codes. Secondly, 
interaction with a companion can be the dominate mass loss mechanism in many post-AGB 
stars, either via tidal interaction or a common envelope. 

Effective temperature. It is not clear what temperature to use. There is no 'natural' temper- 
ature for any physical effect, although the transition occurs around an effective temperature 
of T ~ 5000 K (e.g. Schonberner 1981). Even dust formation can cease at different temper- 
atures, depending on the metallicity of the envelope. Vassihadis & Wood (1994) took the 
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transition to occur when the effective temperature is twice the minimum temperature the 

star can reach on the AGB, but no physical reason is given for that. 

Contraction. Contraction cannot be used because the star starts to contract before it leaves 
the AGB. 

Rapid contraction. The criterion of a rapid contraction, with time or with decreasing en- 
velope mass, captures the essence of the transition, but a quantitative value is not easy to 
define. One can use the logarithmic derivative of the stellar radius with envelope mass 

. d\nR . . 

= . (1) 

cilnMenv 

But S changes monotonically in the relevant temperature (radius) range, and it is not clear 
what value should be used, although 6 = 1 might be a natural choice (Frankowski, A. 2007, 
private communication). Alternatively, one can define the transition to occur when the 
magnitude of the second logarithmic derivative of the stellar radius with envelope mass 

reaches its maximum value. This occurs when the contraction changes from the AGB type 
behavior to the post- AGB one. Examining some models show that this occurs at an effective 
temperature of ~ 7000 — 9000 K, which is hotter than the usually assumed transition point 
(Schonberner & Blocker 1993), e.g., Szczerba et al. (2001) and Tylenda ct al. (2001) who 
listed G and K stars as post AGB stars. The criterion proposed in the next section includes 
in it the beginning of the rapid contraction with envelope mass, and does it with a quantita- 
tive physical definition. In any case, it seems that the criterion of maximum C is similar in 
some aspects to the criterion suggested in the next section (but not the quantitative value 
of transition) . 

Topology in theU — V plane. One defines the quantities V = Anr^p/Mr and U = GMrp/rP, 
where Af,. is the mass inner to radius r in the star, and the other symbols have their usual 
meaning. One can draw the structure of the star in the U — V plane. As discussed in detailed 
by Sugimoto & Fujimoto (2000), a structural curve with a loop corresponds to a giant-like 
structure, while that without a loop corresponds to a dwarf-like structure. Equivalently, one 
can consider the variation oiW = V/W within the star. If it changes (inside the star) mono- 
tonically then there is no loop. We can try to apply this criterion to the post AGB star. For 
a very low envelope mass Mr is constant in the envelope, and W = GM'^ /Ai^r^P ~ {r'^P)~^. 
I examine the structure of the model with the envelope mass of 5.74 x 10~^ from Soker 
(1992), that has an effective temperature of 7200ir. The envelope pressure profile changes 
from P ~ in the range r < 3/?© to P ~ 7--3-7 jj^ ^j^e range 4 < r < 2QRq. Namely, W 
increases with r for r < 3i?0, and decreases with increasing r in the range 4 < r < 207?©. 
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A loop in the U — V plane does exist for this model. Therefore, the loop in the U — V 
plane disappears too late in the post AGB evolution. An interesting property is that V 
and U are related to dynamical and thermal properties of the envelope, however, it is not 
straightforward to connect them to relevant properties of the post- AGB envelope. 
Disappearance of the envelope convective zone. Although the convective zone becomes thin- 
ner and thinner as the star shrinks, it disappear only when the star is very hot (e.g., Soker 
1992). 

End of shell burning. This criterion has not been suggested, but it is listed here for the 
sake of clarity. It is not a good criterion as nuclear burning ceases when the star is hot (e.g. 
Harpaz & Kovetz 1981; Kovetz & Harpaz 1981; Schonberner 1981). 



2. THE PROPOSED CRITERION 

When the envelope mass is low, the Kelvin-Helmholtz time of the envelope is given by 
(Schonberner & Blocker 1993) 

GM, f 47rrV(r) ^ 
TKH-cnv = — ^ / dr, (3) 

where Mc is the core mass, L is the stellar luminosity, R is the stellar radius, and p is the 
density in the envelope. On the upper AGB the density profile in most of the envelope 
(beside regions very close to the core that contain very little mass) can be approximated 
by p oc (e.g., Soker 1992). For the inner radius of the envelope, and in particular for 
the convective part, we can take ro ~ l-R©, and so Ps = ln(-R/ro) ~ 6 in the integration of 
equation I defined a parameter that depends on the exact structure of the envelope. 
For the response of the envelope alone, the internal energy of the envelope should also be 
considered. This will reduce the required time to supply energy or to remove energy, and 
will reduce somewhat the effective value of /3s- I will therefore scale it with /?s = 5. Scaling 
the different variables of upper AGB stars gives 

"""-^ " K^) (^) (^) (h^) (^) 

where Mgnv is the envelope mass. 

For the relevant dynamical time I take (Gpav)^^^^, where pav is the average density of 
the entire star. Scaling with typical number gives 

/ M / R 
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During the evolution along the AGB before the star starts to contract in radius, the 
luminosity increases and the mass decreases, such that Tdyn increases and rxH-cnv decreases. 
Therefore, their ratio Q = Tdyn/rKH-env increases. For evaluating the value of Q on the upper 
AGB, when the envelope mass is low, I take = M in equation (jlj). This gives for the 
upper AGB 

. ^ . 0.1 f f V (JL\ "'^ " (6) 
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The thermal time TKH-env is less than an order of magnitude longer than the dynamical 
time during this late AGB phase. Soker & Harpaz (1999) noted that this relatively short 
thermal time must result in a strong irregular behavior of the envelope because dynamical 
motions, such as pulsations and convective motion, can cause large thermal perturbations in 
the envelope. Therefore, a large value of this parameter, i.e., Q > 0.1, can be related to a 
highly inhomogeneous mass loss process, as well as to a high mass loss rate. Soker & Harpaz 
(1992) argued that the characteristics of AGB stellar pulsations depend on the thermal and 
dynamical time scales. They used the thermal time scale of only the upper envelope, and 
for the dynamical time scale they took the pulsation period. 

The star starts to contract before it leaves the AGB. The initial contraction is slow, and 
the stellar radius during the early contraction phase can be approximated by 

Rc^rJj^X, (7) 

where Rm and Menv-m are the stellar radius and envelope mass when the contraction starts. 
For the model used by Soker (1992) relation ([7]) holds for an envelope mass of 0.001 < Mgnv ^ 
0.1 M0 with (5 ~ 0.2 for most of the time. Then 6 increases more and more rapidly until 
it reaches a very large value when the star contracts by two order of magnitude for a tiny 
change in the envelope mass (Schonberner 1983). Qualitatively similar behavior is found for 
other core masses, but at different envelope masses (Schonberner 1983; Frankowski 2003). 
Using equation ([7]) to express the stellar radius in equation ([6]) gives for the contracting- AGB 
phase 
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(8) 

During the contracting- AGB phase the luminosity and mass do not change much, and 
the derivative of equation ([H]) can be written as 

~ rflnMen, ~ 2 dlnMenv' 



- 6 - 



Along the entire AGB AQ is negative (beside temporal variations, e.g., after thermal pulses). 
It is very positive during the fast contraction along the post-AGB track (again, beside 
temporal variations). The transition from AQ < to AQ > can mark the beginning 
of the post-AGB phase. (Note that the envelope mass decreases with time, and therefore 
when AQ < then Q increases with time.) Namely, the star is said to terminate the AGB 
when Q is at its maximum value; this occurs after the contraction started and Q = Qc- 

If there is no change in the density profile then din (3s/ din M^^^ = and AQc changes 
sign when 6 = 0.4. This is when more or less the rapid contraction starts. However, during 
the contracting- AGB phase the density profile becomes steeper (e.g., Soker 1992), and Ps 
increases slowly, so that din Ps/ din Mem < 0. On the other hand, the envelope convective 
zone, which might be more relevant to many processes influencing the mass loss process, 
becomes concentrated in the outer region and the effective value of jSg might decrease. Over 
all, I suggest to ignore the structural factor Ps, and to mark the transition when 6 = 0.4, 
i.e., when d\n{Rc)/d\nMenv = 0.4 and AQ = 0. 

Alternatively we can use the criterion that 

d In Menv d In Mem 

is equal zero. Namely, XKH-env has its minimum value. This occurs when 5 ~ 1. 

In the model presented in Soker (1992) 6 = 2/5 and 6=1 when the stellar radius is 
R = 150Rq and R = IIORq, and the effective temperature is T = 4400 K and T = 5100 K, 
respectively. The value of C defined in equation reaches its maximum at T ~ 7000 K. In 
a model with a core mass of Mc = 0.67Mq from the solar metallicity track of Vassiliadis & 
Wood (1994), the maxima in Q and Qkh as calculated by Frankowski (Frankowski, A. 2007, 
private communication) are reached when the stellar radius is R = 120Rq and R = 85Rq, 
and at an effective temperature of T = 5200 K and T = 6100 K, respectively. The value of 
C has its maximum at a temperature of T = 8900 K in that model. As mentioned in section 
1, the C-criterion gives the transition at a too high temperature, after mass loss has decline, 
but still captures most aspects of the transition. I prefer the Q-criterion or the QKH-criterion 
because they seem to have more of a physical implication to the mass loss process. 



3. SUMMARY 

I suggest to theoretically define the transition from the AGB to the post-AGB phase, 
namely, the termination of the AGB, by using the envelope thermal time scale (Kelvin- 
Helmholtz time) and the stellar dynamical time. The criterion does not refer to the short 
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time-scale variations occurring on the AGB and post-AGB, e.g., thermal pulses (helium shell 
flashes) and magnetic activity, but refers to the time average of stellar properties. All other 
alternatives for the transition must use average values as well; only a criterion based on 
the envelope mass alone does not need a time average, but no such a criterion has been 
suggested. 

As the star expands along the AGB and loses mass the thermal time scale decreases and 
the dynamical time increase. On the upper AGB the two time scales become comparable. 
This implies that dynamical processes, such as pulsation and convective motion, can influence 
the thermal state of the envelope. This is likely to influence the mass loss process. The 
thermal time scale continues to decrease even after the star starts its contraction. The 
thermal time scale starts to increase at about the same evolutionary point where the rapid 
contraction starts. 

I suggest to deflne the termination of the AGB when Q (deflned in eq. [6]) reaches it 
maximum value. This occurs during the contraction part of the AGB, when Q = Qc (eq. 
E]), and the transition occurs when AQc = (deflned in eq. M)- Alternatively, using only 
the thermal time scale, the transition can be deflned when the thermal time scale starts to 
increase with decreasing envelope mass. Namely, when the value of rxH-env is at its minimum 
and AQkh = (eq. [10]). 

The criterion proposed here has several advantages. 

1. It uses well deflned properties of the envelope. These are simple to derive with a stellar 
model calculated without the inclusion of pulsation or magnetic activity. If these are 
included, time average must be introduced. 

2. The transition can be deflned to occur at a well deflne evolutionary point (when aver- 
aged over short time scales variations), when AQ = (or, alternatively, AQkh = 0). 

3. It is based on properties that are closely related to the mass loss process: The dynamical 
time and the thermal properties of the envelope. However, the theoretical transition 
occurs when Q is at its maximum value. Therefore, the high mass loss rate and 
inhomogeneities are likely to continue into the post-AGB phase for a short time. 

4. It contains previously proposed criteria, (a) Pulsation: via the dynamical time scale, 
(b) Rapid contraction: as the change in the behavior of TKH-env occurs at about the 
same time rapid contraction stars, (c) Mass loss rate: The mass loss process is related to 
the dynamical and thermal time scales, (d) Effective temperature: The value AQ = 
is reached at T ~ 4000 — 6000 K, a temperatures range that was used before to mark 
the zero age pos-AGB phase. 
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Although the criterion proposed here has its own merit in the theoretical study of 
stellar evolution, it goes beyond a pure academic exercise. The criterion suggests that the 
involved time scales have a role in determining the evolution of the envelope, mainly as they 
are important parameters in the mass loss process, e.g., via pulsation properties. Future 
theoretical studies will have to examine in more detail how these time scales affect the 
mass loss process. Also, the suggested theoretical criterion must be examined with detailed 
numerical simulations, to verify that the criterion AQ = 0, or AQkh = 0, do not give wrong 
results. The numerical study should also determine which of these two similar criteria fit 
observational parameters better, hence connecting theoretical studies with observations. 

This research was supported by the Asher Space Research Institute in the Technion. 
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